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Popular Summary 

Profiles of hydrometeor characteristics were estimated using an iterative retrieval 
algorithm. The algorithm minimized the differences between forward calculations and 
observed radar and radiometer observations from the ER-2 aircraft obtained during 
C AMEX-3. The advantages of this retrieval algorithm are the use of high frequency 
channels (89 to 340 GHz) to provide details of the frozen hydrometeors and combining 
radar and radiometer observations. 


The retrieval was performed on anvil, convective, and quasi-stratiform cloud types 
associated with the outer eyewalls of Hurricane Bonnie (August 1998). The retrieval 
results were qualitatively validated using observations from two independent 
measurements observing the same scene. The validation data corroborated the retrieved 
contents and particle size distributions. 

This work shows the importance of including the high frequencies when attempting to 
estimate frozen hydrometeor characteristics using radar and radiometer data. Indeed, 
neglecting the high frequencies yielded acceptable estimates of the liquid profiles, but the 
ice profiles were poorly retrieved. The retrieved profiles contain considerable information 
about the cloud structure and hydrometeor size distributions even though they may not 
have unique solutions. The retrieved profiles are important for improving global change 
models and cloud resolving models. 
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Abstract 

Information about the vertical microphysical cloud structure is useful in many modeling 
and predictive practices. Radiometers and radars are used to observe hydrometeor 
properties. This paper describes an iterative retrieval algorithm that combines the use of 
airborne active and wideband (10 to 340 GHz) passive observations to estimate the 
vertical content and particle size distributions of liquid and frozen hydrometeors. The 
physically-based retrieval algorithm relies on the high frequencies (> 89 GHz) to provide 
details on the frozen hydrometeors. Neglecting the high frequencies yielded acceptable 
estimates of the liquid profiles, but the ice profiles were poorly retrieved. Airborne radar 
and radiometer observations from the third Convection and Moisture Experiment 
(CAMEX-3) were used in the retrieval algorithm as constraints. Nadir profiles were 
estimated for a minute each of flight time (approximately 12.5 km along track) from an 
anvil, convection, and quasi-strati form rain. The complex structure of the frozen 
hydrometeors required the most iterations for convergence for the anvil cloud type. The 
wideband observations were found to more than double the estimated frozen hydrometeor 
content as compared to retrievals using only 90~GHz and below. The convective and 
quasi-stratiform quickly reached convergence (minimized difference between 
observations and calculations using the estimated profiles). A qualitative validation using 
coincident in situ CAMEX-3 observations shows that the retrieved particle size 
distributions are well corroborated with independent measurements. 
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Abstract 


Information about, the vertical microphysical cloud structure is useful in many 
modeling and predictive practices. Radiometers and radars are used to observe hydrometeor 
properties. This paper describes an iterative retrieval algorithm that combines the use 
of airborne active and wideband (10 to 340 GHz) passive observations to estimate the 
vertical content and particle size distributions of liquid and frozen hydrometeors. The 
physically-based retrieval algorithm relies on high frequencies (> 150 GHz) to provide 
details on the frozen hydrometeors. Neglecting the high frequencies yielded acceptable 
estimates of the liquid profiles, but the ice profiles were poorly retrieved. Airborne radar and 
radiometer observations from the third Convection and Moisture Experiment (CAMRX-3) 
were used in the retrieval algorithm as constraints. Nadir profiles were estimated for a 
minute each of flight time (approximately 12.5 km along track) from an anvil, convection, 
and quasi-stratiform rain. The complex structure of the frozen hydrometeors required 
the most iterations for convergence for the anvil cloud type. The wideband observations 
were found to more than double the estimated frozen hydrometeor content as compared 
to retrievals using only 90 GHz and below. The convective and quasi-stratiform quickly 
reached convergence (minimized difference between observations and calculations using the 
estimated profiles). A qualitative validation using coincident in situ C AMEX-3 observations 
shows that the retrieved particle size distributions are well corroborated with independent 
measurements. 
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1. Introduction 


Knowledge of the vertical microphysical cloud structure is important for many 
aspects of meteorology, such as for determination of precipitation rates, latent heating 
profiles, and for forecasting hurricane intensity (Simpson et al. 1996). In addition, 
hydrometeor profiles are used to improve global change models and cloud resolving models. 
Severe storms or intense rain can also affect, earth-satellite communication transmissions. 
For these reasons, accurate estimates of the vertical profile of liquid and frozen hydrometeor 
particle size distributions are vital to the atmospheric research, meteorological, and 
communications communities. In an effort to estimate precipitation profile information 
despite sparsely-situated ground-based sensors, airborne- and satellite- based remote sensing 
instruments have been employed (Kummerow et al. 2000). 

The challenge of using airborne and/or satellite remote sensors is determining 
the appropriate instruments for the parameter of interest. Infrared instruments provide 
temperature and relative humidity profiles in cloud free regions. Lidars remotely sense 
aerosols, clouds, wind speed and direction, and total precipitable water. However, infrared 
and lidar instruments cannot be used to reliably obtain detailed precipitating hydrometeor 
information. A single channel active microwave radar can only provide one of the 2-4 key 
parameters needed to fully characterize the particle size distribution (PSD) for each of 
its range gates, but radar reflectivity information can be contaminated by hydrometeor 
attenuation. A passive multi-frequency microwave radiometer allows probing into the 
different hydrometeor layers of the clouds and the different channels are sensitive to various 
hydrometeor types (e.g., liquid versus frozen). The high frequencies (> 89 GHz) of the 
radiometer are more sensitive to frozen hydrometeors, while the low frequencies are mostly 
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sensitive to liquid hydrometeors. However, radiometers are limited to sensing vertically- 
integrated information about the hydrometeor structure. In addition, the relationships 
between hydrometeor characteristics and the upwelling brightness temperatures are both 
non-linear and non-unique. 

By combining active radar and passive radiometers, the opportunities to estimate 
hydrometeor profiles and cloud characteristics improve (Marzano et al. 1999). In fact the 
Tropical Rainfall Measuring Mission (TRMM) (Kummerowet al. 2000) was the first satellite 
to have both a radar and radiometer designed to measure rainfall. Several radar-radiometer 
retrieval algorithms have been developed for use with the TRMM satellite (e.g., Olson et al. 
1996; Sauvageot 1996; Viltard et al. 2000). Prior to TRMM, most existing remote sensing 
methodologies for estimating cloud structure independently relied on either radiometer or 
radar observations (e.g., Meneghini et al. 1997). 

Associated with TRMM are calibration/validation field campaigns. One such 
field campaign was the third Convection and Moisture Experiment (CAM EX-3) which 
was based in south Florida during August and September of 1998 (Geerts et al. 2000). 
The Texas-Florida Underflights (TEFLUN-B) field campaign combined resources with 
CAM EX-3, with the purpose of underflying the TRMM satellite. Multiple instruments 
located on ground, on low and high altitude aircraft, and on satellites were used to observe 
convective and hurricane systems. Of particular interest for this work are measurements 
from instruments on the high-altitude ER-2 aircraft that provide a single active channel at 
9.6 GHz and 11 brightness temperature channels ranging from 10.7 to 340 GHz. The higher 
frequency channels are extremely useful for determining and constraining the particle size 
distributions of the frozen hydrometeors (Deeter and Evans 2000) and provide a unique 
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aspect to this work in relation to other combined radar-radiometer retrieval algorithms (e.g., 
Marzano et al. 1999). Two other CAMEX-3 instruments provide data for retrieval result 
validation. 

The retrieval algorithm derived herein minimizes the differences between (active 
and passive) observations and forward calculations based on the iteratively estimated 
hydrometeor profiles. The observations are used to constrain the solution. In Section 2, 
the retrieval algorithm will be detailed. The application of the algorithm to the CAMEX-3 
data is described in Section 3, the corresponding retrieval results and validation in Section 
4, followed by a summary in Section 5. 

2. Retrieval Algorithm 

The retrieval algorithm uses both radar and radiometer observations in the 
retrieval process. The algorithm minimizes the error between observations and calculations 
based on the iteratively estimated profile. Hydrometeor contents and drop size distributions 
for rain, cloud water, and frozen hydrometeors are adjusted for each iteration. The flowchart 
of the algorithm is shown in Eig. 1 . The flowchart shows that the radar and radiometer 
calculations are separate, while the error analysis and profile adjustment scheme are 
combined. Each of the components in the flowchart are described in detail in this section. 

The algorithm is initialized by converting nadir-viewed radar reflectivity profiles 
into estimates of hydrometeor content profiles. The fine (37.5 m) resolution of the radar 
range gates is averaged to 0.5 km vertical slabs that extend from 0 to 18 km. The radar- 
to-microphysical profile algorithm is initialized using preset temperature, pressure, relative 
humidity, and cloud water profiles taken from hydrologically appropriate Goddard Cumulus 
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Ensemble (GCE) (Tao and Simpson 1993) profiles. The hydrometeor content, profiles from 
the radar-to-microphysical profile algorithm are partitioned into liquid and frozen particles 
with Marshall- Palmer (Marshall and Palmer 1948) and Sekhon-Srivastava (Sekhon and 
Srivastava 1970) exponential drop size distributions, respectively. While continuity of the 
precipitation flux across the freezing level is not explicitly checked, the masses obtained 
from the radar reflectivities should have smooth transitions from one level to the next. 

If necessary, hydrometeors in the cloud profile are extended one additional height level 
(above, below, or both above and below) the cloud boundaries in order to get non-zero Z 
calculations at all non-zero observed Z heights. The cloud liquid water, rain, and frozen 
hydrometeor contents and drop size distributions are then adjusted at each iterative stage 
of the retrieval algorithm. 

2.1 Brightness Temperature Calculations 

An efficient radiative transfer (RT) model is required to transform the microphysical 
information into upwelling passive microwave brightness temperatures (Tg) that are then 
compared to observations in the iterative retrieval algorithm. The planar-stratified, 
scattering-based RT model used herein was originally developed by Gasiewski (Gasiewski 
1993) and later modified by Skofronick-Jackson and Gasiewski (Skofronick-Jackson and 
Gasiewski 1995) to allow for five (or more) hydrometeor types (e.g., suspended cloud water, 
rain, suspended cloud ice, snow and graupel). Flexibility exists in that the user can input 
the cloud profile and select observation height, viewing angle, frequency (tested from 6 
to 425 GHz) and polarization. The RT model requires as input instrument specifications, 
vertical profiles of temperature, height, relative humidity, and PSD of the hydrometeors in 
the cloud. 
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Typically RT models that handle hydrometeor scattering using a perturbation 
method (Lenoble 1985), as opposed to the adding-doubling method (van de Hulst 1980), as 
the model used herein does, require more iterations to reach the final brightness temperature 
value than the adding-doubling method. An analytical approximation was implemented 
for this retrieval in order to reduce processing time. The RT perturbation technique sums 


successive orders of scattering: 


M 

T b = Z AT b } 


( 1 ) 


*=0 


where ATy is the clear air solution and M is preset, during the perturbation radiative 
transfer modeling. The n to M successive orders of scattering can be written in the form of: 


AT { g ] 


AT’tn+i) A ri n+2) A Tt +1) AT b 

1+ + TTT7T + 


AT, 


(n) 


B 


AT 


B 

,(n+l) 

B 


(M) 


. . . A7 1 , 


(n+l)' 


B 


AT 


(n) 


AT 


,(M — 1) 
B 


AT, 


(«) 


B 


If it can be assumed that AT B + ^ j AT B ^ remains fixed for all k > n and M 
analytical expression for Eq. 2 is: 


A7'L n) - 


1 


1 - W 


where 


W = 


AT, 


("+ 1 ) 


7^( n +*) m 
*_B 

7 . ( 11 ) r p{ n ~ 1 ) 

1 B ~ 1 B 


(") 


AT, 


(n) 


B 


(2) 

oo then an 

( 3 ) 

( 4 ) 


The Tg^ are the brightness temperature values at perturbation k. This approximation can 
cut the number of successive orders of scattering by more than 60%. 


It is assumed that a planar-stratified model is acceptable in this work because 
only high-resolution nadir-viewed observations are used in the retrieval and because the 
observations have fairly small footprints. The calculations assume a nadir-viewed ocean 
surface at 300 K with a surface windspeed of 10 m/s. Similar to Marzano et al. (1999), the 
first 500 m above the ocean surface is not adjusted in the retrieval because of difficulties 
separating surface radar return from the hydrometeors at those levels. Another assumption 
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is that, all particles axe spherical in shape and that frozen hydrometeors are solid ice spheres 
and thus will not need dielectric mixing approximations to compute the absorption and 
scattering coefficients. At high frequencies (> 150 GHz) dielectric mixing theories begin to 
breakdown because the particle inclusion size becomes electrically large with respect to the 
wavelength (Sihvola 1989). 

2.2 Radar Reflectivity Calculations 

The radar reflectivity ( Z ) calculations are based on a model described by Jones 
et al. (1997). This code is also flexible because the user can specify instrument viewing 
angle, height, range gate intervals, etc. The radar model described in Jones et al. 

(1997) relies on GCE data for input and models reflectivities as would be seen from the 
Precipitation Radar on the Tropical Rainfall Measuring Mission ( I RMM) satellite (Simpson 
et al. 1996). For this work, the radar code has been revised to simulate the ER-2 Doppler 
Radar (EDOP) instrument Heymsfield et al. (1996). The model described in Jones et al. 
(1997) automatically generates a melting layer; in this work the retrievals rely solely on 
the hydrometeor parameterization as iteratively generated by the retrieval. Using the size 
distributions, the reflectivity and attenuation can be computed at each range gate. 1 he 
backscattering cross section used in obtaining the reflectivity is computed using Mie theory. 
By integrating the reflectivity over a simulated radar beam that models the observation 
radar specifications, the return power can be calculated and compared to the observed data. 

2.3 Error Analysis 

Once the T fl and Z have been computed for the current estimated profile, the 
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error analysis occurs. There are two convergence criteria. One for the T B and one for the Z: 


max(| T Bobs {i) - T Bcalc {i ) |) < S Tg for all i frequencies (5) 

max(| Z obs (j) - Z ealc (j ) |) < S z for all j range gates. (6) 

The St b and Sz are convergence thresholds and are currently set to 10 K and 2.5 dBZ. Two 
convergence criteria are needed so that the brightness temperature and reflectivity errors 
can be independently analyzed to determine the drop size distribution adjustment needed 
to reduce any errors. 

2.4 Hydrometeor Adjustment Algorithm 

The hydrometeor adjustment algorithm is outlined in Fig. 2. 1 he algorithm first 
checks to see if the radar convergence criterion is met. If not the particle size distributions 
are adjusted at the height h m of the maximum error in the Z profile. The adjustment factor 
Y is related to the ratio {Z obs {h m ) / Z ea i c {h m )). The Z-based adjustments are 

A = A old /Y and N 0 = N° ld /Y 3 (7) 

in the exponential PSD equation 

A 1(D) = N 0 e~ AD (mm -1 m -3 ) (8) 

where £>, A -1 are in mm and N 0 is in mm -1 m -3 . This application of Y preserves the T B 
values so that they do not change much for Z-based adjustments. The starting No and A 
are provided during the initialization procedure previously described. 

When the convergence criterion for the T B is not met, a similar adjustment is made 
to reduce the differences in the observed and calculated T B values. Frequencies where the 
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temperature difference exceeds 10 K are used to determine where to adjust the hydrometeor 
profile. Each frequency has a weighting function that peaks at a height dependent on 
temperature and hydrometeor characteristics. The hydrometeor profile is adjusted at the 
height(s) and ±500 m of where the weighting function(s) peaks for the frequency(ies) with 
error exceeding 10 K. If the radiometer frequencies are less than 89 GHz, the rain and 
cloud water drop sizes are adjusted; if above 89 GHz the frozen hydrometeor particle sizes 
are adjusted. For the 183.3 GHz channels, only the 183. 3±1 GHz channel is used to adjust 
rain and cloud water drop sizes. Errors in the 89, 183. 3±3 and 183. 3±7 GHz channels are 
usually corrected through PSD adjustments due to errors in other channels. 

If the weighting function peak does not fall within the cloud boundaries, the 
particle size adjustments are made to the lower, middle, and upper third of the cloud 
volume based on the frequency. For 37 and 150 GHz, particle sizes within the lower third 
of the cloud liquid and frozen hydrometeor layers (respectively) are adjusted. For 183. 3±1 
and 220 GHz, the particle sizes within the middle third of the cloud liquid and frozen 
hydrometeors (respectively) are adjusted. Finally for 340 GHz, the particle sizes in the 
upper third of the frozen hydrometeor layers in the cloud are adjusted. 

The adjustment parameters are a function of the ratio (Tsobs(k) /TB ca ic{k)) denoted 
by X where k is the index of the frequency at which an error occurs. The 7’s-based 
adjustments are 

A = A old /X and N 0 = N° ld / X 7 (9) 

in the exponential particle size distribution (PSD) equation (8). This application of X 
preserves the Z values because Z is proportional to D 7 No (Z cc J D e N(D)dD). 

The adjustments to the size distributions of a specific height level and hydrometeor 
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type are only changed once per iteration. This eliminates the scenario where one adjustment 
for a frequency or range gate would cancel out the adjustment for a different frequency. 
Also the adjustments do not preserve the (bulk) density of the hydrometeors because 
the bulk density M is proportional to D 4 N 0 . While we have been able to meet our 
convergence criteria for anvil, convective, and stratiform cloud cases, it is likely that further 
improvements to the adjustment algorithm will allow us to tighten the convergence criteria 
thresholds. Suggested improvements will be discussed in the Summary section. 

3. Application to CAMEX-3 Data 

Observations from CAMEX-3/TEFLUN-B (Geerts et al. 2000) are used in the 
retrieval algorithm. While CAMEX-3 had broad-based instrumentation on multiple 
platforms including aircraft, ground, ship, and satellite, this work focuses on the data 
from four instruments on the ER-2 aircraft and one instrument on the DC-8 aircraft. On 
board the ER-2, flying at an altitude of 20 km, the instruments of interest for this work 
are the Millimeter-wave Imaging Radiometer (MIR) Racette et al. (1996), the Advanced 
Microwave Precipitation Radiometer (AMPR) Spencer et al. (1994), the ER-2 Doppler 
Radar (EDOP) Heymsfield et al. (1996) and the National Polar-orbiting Operational 
Environmental Satellite System (NPOESS) Aircraft Sounder Testbed - Microwave (NAST- 
M) (Blackwell et al. 2001). These ER-2 instruments measure atmospheric hydrometeors 
in the microwave region of the electromagnetic spectrum. The MIR observes at 89, 150, 
183.3±1, ±3, ±7, 220, and 340 GHz, while the AMPR observes at 10.7, 19.35, 37, and 
85.5 GHz. The EDOP is an active radar sampling at 9.6 GHz with a range gate interval 
of 37 meters. The NAST-M, used only for validation purposes, is a passive microwave 
spectrometer with 17 channels near the oxygen absorption lines at 50-57 GHz and 
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118.75 GHz. The NAST-M channels are independent of any channels used in the retrieval. 
On the DC-8 aircraft flying at 12 km altitude, there is the Particle Measuring System, Inc. 
(PMS) 2DC probe as part of the Cloud and Aerosol Particle Characterization (CAPAC) 
suite of instruments. The PMS 2DC probe provides in situ observations of particle size 
distributions used to validate the retrieved size distributions. 

For this work, only the nadir or near-nadir signatures are used. The data sets 
from the MIR, AMPR, and EDOP have been analyzed and colocated. In an effort to match 
each radiometer channel and the corresponding EDOP radar samples to similar volumes 
of the storm, all data values were simulated to match the lowest resolution channel of all 
the instruments. The largest footprint of 2.8 km occurs for the AMPR 10 and 19 GHz 
channels at nadir when the ER-2 is flying at 20 km altitude. Because there is not an 
AMPR pixel directly at nadir, for the 10 and 19 GHz channels the two pixels adjacent to 
nadir were averaged in an effort to simulate a nadir value. In order to transform the MIR, 
EDOP and the higher frequency AMPR data to the lowest resolution, a two-dimensional 
Gaussian-weighted mean for each set of data values was calculated where the weights are 
1.0 at the center and 0.5 at 1.4 km away from the center. Additionally, since the EDOP 
radar does not sample off nadir, it probably misses some of the volume sample that the 
radiometers observe. The EDOP has a 3° beam width so it is seeing about 0.6 km off nadir 
in both directions at the surface instead of the 1.4 km off nadir view of the AMPR 10 and 
19 GHz channels. 

While checking the quality of the match between the instruments an offset between 
the AMPR 85 GHz channel and the MIR 89 GHz channel was observed. Even though the 
quantitative values of the brightness temperatures of these two channels are not expected to 
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match exactly, it is reasonable to assume that the time series of the brightness temperature 
trend for the two channels should be similar. It was determined that the value of the 
offset in time was dependent upon ER-2 altitude and airspeed; therefore, a pointing angle 
difference was hypothesized as the cause of the error. This pointing angle difference was 
consistently determined as 3.3 degrees when using three separate sets of field campaign 
data for MIR and AMPR from the 1998 (CAMEX-3) and 1999 (TRMM Large Scale 
Biosphere- Atmosphere Experiment in Amazonia (LBA)) time period. The cause of this 
error is likely related to slight differences in the orientation of instruments within the ER-2 
aircraft. The correlation of the time series of the EDOP and AMPR data was determined 
to be acceptable so no correction was applied to these data; MIR data was corrected using 
the pointing angle difference along with the ER-2 altitude and airspeed values. 

Figure 3 shows the colocated observed data for Hurricane Bonnie on 26 August 
1998, including EDOP data (upper panel), MIR brightness temperatures (center panel), 
and AMPR brightness temperatures (lower panel). For this image, the ER-2 is flying west 
to east at approximately 32.8° N latitude. The left side of the image corresponds to the 
outer edge of Hurricane Bonnie and the plane flies toward the eye, which is on the right 
side of the image, crossing several rainbands. There are at least three distinct cloud types 
in Figure 3: an anvil cloud region to the left of the image (14:04:00 UTC to about 14:07:30 
UTC) with a high altitude EDOP reflectivity profile and low brightness temperature values 
for the MIR frequencies (indicating cooling from ice scattering). There is a weak outer 
eyewall embedded convective core from about 14:11:00 UTC to about 14:13:30 UTC and a 
quasi-stratiform cloud region from about 14:15:20 UTC to about 14:17:00 UTC. 

The retrieval algorithm focuses on small regions within these three storm types. 
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The selected times span from 14:05:04 — 14:06:05 UTC, 14:12:07 — 14:13:41 UTC, and 
14:15:40 — 14:16:38 UTC for the anvil, convective and quasi-stratiform regions respectively. 
The temporal sampling of the colocated data set is about 3 seconds and has a footprint size 
of 2.8 km at the ER-2 altitude of 20 km. However, since the observations are oversampled, 
the retrieved profiles are between 12.5 and 18.7 km along-track of the flight line, for a total 
of ^ 44 km. An analysis of the complete data set as shown in Fig. 3 is computationally 
prohibitive at this time and would require additional adjustments to the size distribution 
adjustment algorithm for the transitioning cloud type regions. 

4. Retrieval Results 

The iterative retrieval algorithm was used to estimate the vertical hydrometeor 
profiles of the three storm types found within Hurricane Bonnie’s outer rainbands. The 
retrieval performed well, producing extremely good matches between the observed Z 
and calculated Z (see Fig. 4ab). The brightness temperature matches were less than 
the maximum threshold of 10 K for nearly all of the cases (see Fig. 5). The average 
number of iterations required for estimating one profile was 36 with the anvil retrievals 
requiring more iterations than the average number to converge and the quasi-stratiform 
case requiring fewer than the average number of iterations to converge. In general it was 
easiest to reach convergence for the radar reflectivity stopping criterion and more difficult 
for the high frequencies of the radiometer. This is reasonable because the initial estimated 
profiles were obtained by transforming the EDOP reflectivities to hydrometeor content 
profiles. Additionally, the high frequency radiometer channels are quite sensitive to the 
frozen hydromet.eor size distributions, which are not well known, and therefore make it 
difficult to ascertain the proper sizes to obtain a good match between the high frequency 
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observations and calculations. Finally, it must be stated that some of the retrieved profiles 
are not unique solutions for the observed data. In particular, the anvil cases had 2 or more 
slightly different PSD profiles which still met the convergence criteria (see Fig. 6a). It 
is assumed that if the convergence criteria were even more restrictive that the differences 
in the retrieved profiles would decrease. Even though the retrievals may not be unique, 
they still provide good approximations to the actual cloud structure and hydrometeor size 
distributions. 

4-1 Anvil Results 

For the anvil case (14:05:04 — 14:06:05 UTC), output from the estimated profiles 
are shown in Figs. 4c-f (left hand side of each image). Figs. 4c, e show the liquid and frozen 
hydrometeor contents while Figs. 4d, f show the liquid and frozen (melted) hydrometeor 
median diameter (D 0 = 3.67/ A from Eq. 8). The algorithm also produces the number 
density (iVo from Eq. 8) but this value is not plotted because it can be extracted using the 
content M and Eq. 8 in 

J r OO 

( N(D)--pdD (10) 

0 6 

where p = 1.0 g cm -3 for rain and 0.913 g cm -3 for the frozen hydrometeors. 

From Fig. 4e, notice that the frozen hydrometeor content is large across the whole 
anvil region. The variability can be attributed to the difficulties in reaching convergence 
for the anvil region profiles. In the anvil region, the median diameter of the liquid and 
frozen (melted diameter) particles does not exceed 0.16 mm and 0.56 mm, respectively. 
The average (over all heights and UTC times in the anvil cloud) melted median diameter 
is 0.077 mm. The radar and radiometer calculations extend the exponential drop size 
distribution to approximately 2.7 times these median diameters. 
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Fig. 5 shows the differences in calculated and observed 7#. As can be seen from 
Fig. 5 (left hand side), the differences for the 37 and 340 GHz brightness temperatures 
bound the convergence criterion space. The average number of iterations required for 
estimating an anvil profile is 78. The required number of iterations is higher than the other 
cases because of the difficulty in determining appropriate characteristics for the frozen 
hydrometeors. In addition, the assumption of solid ice spheres here rather than specific 
crystalline shapes or fluffy snow may increase convergence time. 

In the anvil region, there is a negative bias where the calculations are too cool 
for the low frequencies and a positive bias for the highest frequencies. It is possible that 
the contents of the monodisperse cloud water and cloud ice initialized by the GCE model 
are not properly estimated. (Note that cloud water content, but not cloud ice content, is 
adjusted in the retrieval.) Quick tests using cloud water and cloud ice from other cloud 
models do not improve the retrieval results significantly. Modifying the retrieval algorithm 
to adjust the cloud water and cloud ice presents several difficulties: (1) the radar data is 
not sensitive to cloud water and cloud ice, (2) cloud water and cloud ice may extend beyond 
the boundaries of the radar returns strong enough to indicate the presence of clouds, and 
(3) how should cloud water and cloud ice be modeled: monodisperse or polydisperse. The 
negative/positive bias may also be caused by inadequate modeling of the ocean surface 
and/or lower altitudes. The calculated high frequency radiometer data is relatively warm 
with respect to the observations because if additional ice is introduced it over cools the 
37 GHz channel. 

Figure 6a shows that the profile retrievals may not be unique. This figure shows 
the retrieved profile for 14:05:04 UTC for three different size distribution adjustment 
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algorithms (modifications to Fig. 2). Algorithm 1 is used for the reported results in this 
paper. All three algorithms met the convergence criteria; however, algorithm 2 tended to 
put more ice at upper altitudes and less at lower altitudes. Although it is not shown, the 
retrieved profiles for the three algorithms are more similar for thicker anvil, convective, and 
stratiform clouds. Although the profiles are different, the columnar ice masses among the 
three algorithms are within 100.0 g/m 2 of each other. 

For anvil clouds the necessity of using the high frequency channels to determine 
frozen hydrometeor characteristics is shown in Fig. 6b. This figure shows what happens 
when the retrieval algorithm does not utilize (i.e. correct for errors in) the high frequency 
channels. For the single profile at 14:05:07 UTC, the solid line in Fig. 6b shows the results 
using all channels and the reflectivities. For this figure, the solid line is considered the 
truth. The dotted line shows the retrievals using only frequencies less than and equal to 
150 GHz. The dashed line shows the results when using only frequencies less than or equal 
to 89 GHz. For the low frequency case (dashed line) the calculated Tg for the 340 GHz 
channels are 60-70 Kelvin warmer than the observations. Clearly much information about, 
frozen hydrometeor content cannot be derived when only the low frequency channels are 
used. 

4-2 Convective Results 

The convective profiles (14:12:07 — 14:13:41 UTC) on the average required fewer 
iterations (25). This is due to the underlying rain profile stabilizing the low frequency 
radiometer data and reducing the impact of the ocean surface conditions. Figures 4 and 5 
(center region) show the retrieved content, median diameter, and 7'g differences. Figure 4 
indicates liquid and frozen contents similar to what would be expected for the EDOP profile 
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(Fig. 4a). The maximum median drop diameter for rain (over all non-zero retrieved rain 
pixels) is 0.58 mm; for the frozen hydrometeors the maximum melted median drop diameter 
is 1.1 mm. It is encouraging that the frozen hydrometeors (Fig. 4e, f center convective 
region) follow the pattern seen in the reflectivity observations. 

The Tb differences in Fig. 5 show some interesting trends. First, the left hand 
side of the convective region shows a similar bias as seen within the anvil region where the 
low frequency calculations underestimate and the high frequency calculations overestimate 
the actual Tb values. This issue is further compounded because the retrieval algorithm 
does not adjust the hydrometeors in the region between the Earth’s surface and 500 m up 
because it is difficult to separate surface radar return from hydrometeors at those levels. 
Second, the 19 and 37 GHz calculations become warmer than the observations to the left 
of the convective region where the heavy rain region is dissipating somewhat. At the same 
time and in the same region, the high frequency Tb differences are converging toward zero. 
This change of Tb differences is likely due to decreasing hydrometeor complexity from left 
to right in the convective region. The decreasing complexity may mean that the solid ice 
spherical particles are adequate and this causes the high frequency Tg differences to be 
reduced on the left hand side of the image. Notice that at 14:12:38 UTC the 220 GHz 
difference did not converge to less than 10 K, instead it converged to 10.4 K. In addition 
the 183.3±1.0 GHz channel differences are 10.2 K at 14:13:20 and 14:13:23 UTC. 

f.S Quasi- Stratiform, Results 

The quasi-stratiform results are quite surprising considering that no melting layer 
is modeled in the RT calculations or estimated profiles and yet the retrieval reached the 
convergence criteria in only an average of 9 iterations. There are two reasons for the quick 
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convergence: (1) a uniform cloud structure and (2) the initialization procedure did a good 
job estimating liquid and the few frozen hydrometeor levels of this storm, indicating that 
the Marshall and Palmer (1948) PSD is reasonably appropriate for this quasi-stratiform 
cloud. Once again content and median diameter profiles are shown in Fig. 4 (right hand 
side) and 7g differences are shown in Fig. 5. The maximum median diameter is 0.39 mm 
and 1.03 mm for liquid and frozen hydrometeors (melted) respectively. 

The frozen hydrometeor content of the quasi-stratiform region shown in Fig. 4 
follows the EDOP reflectivity profile (e.g., more ice on the left hand side of the quasi- 
stratiform region and a peak diameter band at the melting layer height). Note that the 
bright band height increases toward the right as we move toward the hurricane eye. 1 his 
is typical since the temperatures are warmer in the eye. The 7'g differences are revealing: 
there is a trend of improved comparisons between observed and calculated values near the 
center of the quasi-stratiform time segment (14:15:40 — 14:16:38 UTC). At the fringes of 
the time segment there is a steep decrease of rain and/or frozen hydrometeors as shown 
in the EDOP image. Within these highly variable regions it is more difficult to match the 
observed and calculated 7g. 

4-4 Validation of Retrieval Results 

Since the true microphysical profile cannot be measured at the spatial and 
temporal scales appropriate for the nadir-viewed ER-2 aircraft brightness temperature 
observations, qualitative verification of the retrieved profiles will have to suffice. Two 
different sources of CAMEX-3 coincident information are used to validate the results: 

(1) the passive brightness temperature oxygen band channels on the NAST-M (Blackwell 
et al. 2001) and (2) in situ observations at ~ 12 km from the DC-8 aircraft CAPAC 2DC 
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PMS probe. 


First the 50.3±0.09 and 118.75T3.5 GHz NAST-M channels were compared to 
calculations at those same frequencies using the retrieved profiles. These two channels 
represent the NAST-M channels furthest away from the oxygen line centers and thus more 
sensitive to the hydrometeors in addition to the oxygen. The NAST-M instrument is 
onboard the ER-2 sampling the same (nadir- viewed) scene as the EDOP, MIR, AM PR. Since 
the NAST-M channels were not resampled to the grid of the combined EDOP-MIR-AMPR 
observations, there were fewer overlapping samples. Each of the calculated Tb values using 
the retrieved profile for these overlapping UTC times was less than 10 K different from 
the NAST-M values. This means that the retrieved profiles are consistent (in terms of the 
convergence criteria) with I'b values of the outermost wings of these oxygen bands. 

The second validating data set is from the CAPAC suite of in situ microphysical 
measurement instruments observing from the DC-8 aircraft during CAMEX-3; the PMS 
2DC probe was the primary instrument of interest here. The DC-8 flies at a height of 
about 12 km (8 km below the ER-2) and for this UTC segment was well correlated with 
the ER-2 flight path. Since 12 km is above the cloud tops of the convective and stratiform 
regions retrieved herein, the 2DC microphysics data can only be used to validate anvil ice 
characteristics. The 2DC probe images particles within its field of view and then processing 
produces the number of particles in specified bin sizes within a centimeter-cubed volume of 
air. To compare these measurements with the retrieved exponential PSDs, the number of 
particles in the centimeter-cubed volume were divided by the interval size of the diameter 
bin sizes. This provides a numbers considered equivalent to N(D) for each of the diameter 
bin sizes. The retrieved No were multiplied by e -1 to provide a value of N(D) at D = Do- 
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Only the retrieved values near the 12 km height level were used. Figure 7 shows the 
comparison of PMS 2DC data for 14:00 to 14:06 UTC (the DC-8 was underflying the ER-2 
for these UTC times) to the retrieved anvil ice at 11.5 t.o 12.5 km. This figure shows that 
the retrieved values (symbols for each height) are very similar to the PMS 2DC data (lines 
for each UTC time). This qualitative validation indicates that the anvil ice retrievals are 
close to the in situ observations. 

5. Summary 

This paper has provided estimates of precipitation profiles and frozen hydrometeor 
profiles when retrieving using wideband radiometer observations plus radar observations. 
Profiles of hydrometeor characteristics were estimated using an iterative retrieval algorithm. 
The algorithm minimized the differences between forward calculations and observed radar 
and radiometer observations from the ER-2 aircraft obtained during CAMEX-3. The 
advantages of this retrieval algorithm are (1) the use of high frequency channels to provide 
details of the frozen hydrometeors and (2) combining radar and radiometer observations. 

Contents and particle size distributions for spherical rain, cloud water, and frozen 
hydrometeors were estimated for profiles extending to 18 km with vertical spacing of 0.5km. 
The retrieval was performed on anvil, convective, and quasi-stratiform cloud types. The 
anvil cloud type required the most iterations in order to resolve the unknowns related to the 
characteristics of the frozen hydrometeors. The quasi-stratiform region met the convergence 
criterion the quickest because the initialization procedure used drop sizes more applicable 
to stratiform cloud types. 

The retrieval results were qualitatively validated using observations from the 
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NAST-M on the ER-2 and the CAPAC 2DC probe on the DC-8 aircraft. The brightness 
temperatures of the outermost wings of the 50-60, and 118 GHz oxygen bands as measured 
by the NAST-M, were within the convergence criteria (less than 10 Kelvin from the 
calculated brightness temperature values). Likewise, observations from the PMS 2DC probe 
on the DC-8 aircraft flying at 12 km above the Earth’s surface were used to validate the 
particle size distributions of the anvil region retrievals. The retrieved anvil particle sizes 
and number densities matched the measured 2DC probe microphysics well. Unfortunately 
the DC-8 was above the cloud tops for the convective and quasi-stratiform regions and no 
validation using the PMS 2DC data could be performed for them. 

The retrieved profiles contain considerable information about the cloud structure 
and hydrometeor size distributions even though they may not have unique solutions. The 
retrieved profiles are important for improving global change models and cloud resolving 
models. Furthermore this work shows that high frequency microwave channels (> 150 GHz) 
information needed in order to define the frozen hydrometeor characteristics found at the 
upper altitude levels of a cloud. 
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Figure 1: The retrieval algorithm flowchart. 


Figure 2: The size distribution adjustment flowchart. 

Figure 3: The EDOP radar reflectivity (upper panel), the MIR brightness temperatures 
(center two panels) and the AMPR brightness temperatures (lower panel). 

Figure 4: The (a) EDOP radar reflectivity, (b) calculated reflectivity for the retrieved 
profile, (c) retrieved liquid water content, (d) retrieved rain median diameter, (e) retrieved 
frozen content, and (e) retrieved frozen hydrometeor median melted diameter. 

Figure 5: The calculated minus observed brightness temperature differences for low 
frequencies (upper panel), high frequencies (middle panel), and 183.3 GHz channels (lower 
panel). 

Figure 6: Single profile retrieval results for (a) different hydrometeor adjustment retrieval 
algorithms and (b) for a reduced microwave channel set. 

Figure 7: CAPAC PMS 2DC in situ observations from the DC-8 at 12 km altitude for 14:00 
— 14:06 UTC (lines) and retrieved anvil N(D ) for D = Do at 11.5 km (diamonds), 12.0 km 
(triangles), and 12.5 km (squares). 
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Figure 1 : The retrieval algorithm flowchart. 
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Gail Skofronick-Jackson, et al. 

Page 27 



MIR, EDOP, and AMPR - August 26, 1998 


14:04:02 


14:07:46 


14:11:30 
Time (UTC) 


14:15:14 


60.0 

M 

50.0 

CD 

40.0 

“O 

30.0 

>N 

20.0 

> 

10.0 

o 

0.0 

0) 

-10.0« 


O' 



+89 GHz 
*150 GHz 
<>220 GHz 
A 340 GHz 


+183±1 GHz 
*183±3 GHz 
a 183±7 GHz 


14:18:57 


Figure 3: The EDOP radar reflectivity (upper panel), the MIR 
brightness temperatures (center two panels) and the AMPR brightness temperatures(lower panel). 
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Figure 4: The (a) EDOP radar reflectivity, (b) calculated reflectivity for the retrieved profile, (c) 
retrieved liquid water content, (d) retrieved rain median diameter, (e) retrieved frozen content, and 
(e) retrieved frozen hydrometeor median melted diameter. 
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Figure 5: The calculated minus observed brightness temperature differences for low frequencies 
(upper panel), high frequencies (middle panel), and 183.3 GHz channels (lower panel). 
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Figure 6: Single profile results (a) for three different hydrometeor adjustment retrieval algorithms 

and (b) for three microwave channel sets. 
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Figure 7: CAPAC in situ observations from the DC-8 at 12 km altitude for 14:00 to 14:06 UTC 
(lines) and retrieved anvil N(D) for D = D 0 at 1 1.5 km (diamond shapes), 12.0 km (triangles), and 

12.5 km ( squares). 
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